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ABSTRACT 

The new 7-ray burst (GRB) mission Swift has obtained pointed observations of several classical 
novae in outburst. We analyzed all the observations of classical novae from the Swift archive up to 
30 June, 2006. We analyzed usable observations of 12 classical novae and found 4 non-detections, 3 
weak sources and 5 strong sources. This includes detections of 2 novae exhibiting spectra resembling 
those of Super Soft X-ray binary Source spectra (SSS) implying ongoing nuclear burning on the white 
dwarf surface. With these new Swift data, we add to the growing statistics of the X-ray duration and 
characteristics of classical novae. 

Subject headings: stars: individual (V574 Pup, V382 Nor, V1663 Aql, V5116 Sgr, V1047 Cen, V476 
Set, V477 Set, LMC 2005, V1494 Aql, V4743 Sgr, V1188 Sco) — stars: novae - 
X-rays: stars 



1. INTRODUCTION 

The explosions of Classical Novae (CNe) occur in close 
binary star systems where mass is accreted onto a white 
dwarf (WD) from a low-mass secondary star. The ac- 
creted material gradually becomes degenerate, and when 
temperatures become high enough and the pressure at 
the bottom of the accreted envelope exceeds a certain 
value, a thermonuclear runaway results. Enough energy 
is deposited in the accreted material to eject some frac- 
tion from the WD. The outburst can last several months 
to years, and CNe become active in X-rays at some 
time during their outburst (e.g., iPietsch et al.l 120051 : 
iPietsch et al.l l2006) . The relatively slow evolution makes 
CNe ideal targets for Swift, as no accurate times of ob- 
servation need to be scheduled and they can be observed 
with great flexibility. 

The evolution of X-ray emission starts soon after out- 
burst when the expanding nova envelope is still small and 
dense. The energy pe ak during this earl ie st time, called 
the "fireball" phase (jShore et alJ Il994t ISchwarz et all 
l200lh . is expected to occur at X-ray wavelengths. This 
phase is extremely difficult to observe since it is pre- 
dicted to last only a few hours after the beginning of the 
outburst so that it is over before the nova can be dis- 
covered in the optical. As the envelope increases in size 
and cools, the opacity increases, shifting the radiation 
towards wavelengths longer than the hydrogen absorp- 
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tion edge at 13.6eV (911 A). The expanding, cooling en- 
velope gives rise to the large and rapid increase in the 
visual luminosity of the nova. Around the time of visual 
maximum, X-rays from the underlying WD are trapped 
inside the large column density of the ejecta. However, 
shocks may develop within the early expanding ejecta 
producing a hard, sub-Eddington lumi nosity, X-ray spec- 
tral energy distr i butio n (e.g., V382 Vel. lOrio et alJ200ial 
iMukai fc Ishidal 120011 ). As the shell continues to ex- 
pand, the density and opacity drop, and X-rays emitted 
from the surface of t he WD eventually become visible 
(|Krautter et al.lll996h . Soft X -ray emission originates 
from hot layers produced by the nuclear burning of the 
remaining accreted material on the WD. This material 
burns in equil ibrium at a high but constan t luminosity, 

« LEddingt on 

(|Gallagher fc Starrfield! fl978l) . The spec- 
tral energy distribution resembles t hat of the super-soft 
X-ray sources (S SS), (e.g. Cal 83. IPaerels et al.l I200H : 
lLanz et al.ll2005h . 

Recent observations of novae during this phase show 
that the flux can al so be highly va ri able on short 
times cales (see, e.g. iNess et ail 120031 : lOsborne et al.1 
2006). The duration of this SSS phase is pr edicted to 
be in versely proportional to the WD mass (|Starrfieldl 
119921). but a revie w of the ROSAT X-ray sky survey 
(jOrio et al.l l2001bh showed that the timescale is much 
shorter than predicted, implying either that the masses of 
WDs in CNe are much higher than commonly assumed, 
or that the X-ray turn-off is a function of more than 
the WD mass . A d ifferent approach was proposed by 
iGreiner et alJ (|2003l ) who found from the available data 
that systems with shorter orbital periods display long 
durations of supersoft X-ray phases, while long-period 
systems show very short or no SSS phase at all. They 
speculate that shorter periods may be related to a higher 
mass transfer rate (e.g., by increased irradiation) and 
thus to the amount of material accreted before the ex- 
plosion. 

Another source of X-ray emission from novae are 
emission lines from material that has been radiatively 
ionized. These lines are expected to be present during 
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the entire evolution, but are usually only observed once 
the nuclear fuel is consumed and the bright continu- 
ous SSS spectrum (that outshines the emission lines) 
fades. The ejecta quickly recombine until a collisional 
equilibrium is reached. This equilibrium reflects the 
kinetic temperature distribution of th e plasma and th e 
elemental abundances can be derived (|Ness et al.ll2005l) . 
This phase is called t he nebular phase, and is also seen 
at other wavelengths (| Shore ct al. 2003). In some cases, 
X-ray emission from lines has been observed earlier, 
either prior to the SSS phase, superimposed on the SSS 
spectrum (e.g., V1494Aql), or during times of extreme 
variability when the SSS emis sion suddenly d eclines 
before recovery (e.g., V4743 Sgr. fNess et al.ll2003t) . 

X-ray observations provide important clues to the 
properties and dynamics of novae. The evolution of the 
SSS depends on the WD mass, the mass loss rate (via 
radiatively driven winds or common envelope mass loss), 
the amount of ejected material, the amount of material 
remaining on the WD, the binary separation, and the ex- 
pansion velocity of the ejecta. Analysis of the X-ray emis- 
sion can also provide insight into the composition of the 
ejecta. Unfortunately, only a few novae have been stud- 
ied extensively in X-rays with much of the recent infor- 
mation coming from Chandra and XMM-Newton. Now, 
with Swift, we have the possibility of studying novae in 
large numbers in order to assess statistical properties. 
Swift also allows for coordination of higher- resolution ob- 
servations with Chandra and XMM-Newton since it has 
always been difficult to predict the brightness in X-rays 
and to estimate the optimal exposure time for X-ray ob- 
servations with the gratings. 

In this paper we discuss the Swift XRT instrument 
briefly in Sect. 12.11 Since interstellar absorption plays 
an important role in Galactic novae, we discuss the ex- 
pected effects on the detection of SSS in Sect. 12.21 In 
Sect. 131 we sort the observations by non-detections, de- 
tections of weak sources, and strong sources based on our 
extraction statistics detailed out in the appendix section 
IA.2I We discuss our results in Sect. |4] and expand on the 
type of X-ray emission that was detected in the Swift ob- 
servations. Finally, Appendix IA.1I provides background 
information for each nova. 

2. Swift OBSERVATIONS 

2.1. The instrument 

Our entire dataset has been o btained with the X - 
ray Telescope (XRT) aboard Swift (|Burrows et alj|2005l ). 
The XRT instrument is a CCD detector behind a Wolter 
Type I grazing incidence mirror consisting of 12 nested 
shells. The field of view covers 23.6' x 23.6', imaged on 
a detector with 600 x 600 pixels, thus each pixel cor- 
responds to 2.36". The point spread function (PSF) 
can be parameterized, and source radii of 10 pixels and 
5 pixels include 80.5 per cent and 60.5 percen t of the to- 
tal energy, respectively (jMoretti et al.|[2004T ). The posi- 
tional accuracy is 2.5", or one pixel. The energy range 
covers 0.2 - lOkeV with the Full Width at Half Maxi- 
mum (FWHM) energy resolution varying from ~ 50 eV 
at 0.1 keV to ~ 190 eV at lOkeV. In Fig. Q] we show 
a comparison of effective areas of Swift, Chandra, and 
XMM-Newton. In general, the Swift XRT has a smaller 
effective area than the Chandra/ ACIS-S except at low 
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Fig. 1. — Comparison of XRT effective areas of Swift with other 
X-ray instruments that are sensitive in the same spectral region. 
None of the Chandra or XMM-Newton gratings has a higher effec- 
tive area than the XRT. The CCD detectors: ACIS-S (Chandra) 
and EPIC (PN/MOS1; XMM-Newton) have higher effective areas. 
The spectral resolution of the XRT is similar to either the ACIS or 
EPIC detectors, but the spectral resolution of the gratings is much 
superior. 



2.0xio 22 




2x10 5 3x10= 4x10 5 5x 
Effective temperature (K) 



6x10 5 



Fig. 2. — Contour plot of the logio Swift counts predicted with a 
1-ksec observation from a SSS a a function of temperature (Cloudy 
model) and interstellar extinction (Njj) for a source at a distance 
of f kpc. 

energies (long wavelengths) where the Chandra ACIS- 
S suffers from large calibration uncertainties. The EPIC 
detectors PN and MOSl (XMM-Newton) have 50 percent 
and 20 percent higher effective areas, respectively. The 
XRT, ACIS, and EPIC instruments have similar spec- 
tral resolution. The XRT effective area is larger than 
that of the gratings aboard Chandra (LETG and MEG) 
and XMM-Newton (RGS1 and RGS2), but the XRT 
spectral resolution is much lower. Fig. [T] demonstrates 
that the XRT is an ideal instrument for exploring the 
emission level of novae using reasonable exposure times 
(~ 3 — 6ksec). Once a nova is detected with Swift and 
found sufficiently bright, additional observations can be 
requested with the high-resolution grating instruments 
aboard Chandra and XMM-Newton to obtain detailed 
spectral information. The detector was operated in Pho- 
ton Counting (PC) mode which provides 2D imaging, 
spectral information, and 2.5-sec time resolution. 

2.2. Observability of the Super Soft Source phase in 
Classical Novae 
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Unfortunately, many Galactic CNe have large extinc- 
tion rates which has a pronounced effect on the observ- 
ability of the SSS emission. To il lustrate this we at- 
tenuate a series of Cloudy models ([Ferland et al.|[l998f) 
through different amounts of H I column densities and 
determine the resulting S wift counts. T he models use 
non-LTE atmospheres by iRauchl |l997) for planetary 
nebula nuclei with effective temperatures in the range of 
(1 — 6) x 10 5 K (8-50eV) as the input source. The model 
luminosity was fixed at 10 38 erg sec -1 , which is approxi- 
mately the Eddington luminosity for a 1-M© star. The 
spectral energy distributions were attenuated by H I col- 
umn densities ranging from (f — 20) x 10 21 cm -2 and then 
converted to Swift net count rates using the effective area 
of the XRT and scaled to a 1 ksec exposure at 1 kpc. The 
logarithms of the expected total XRT counts are shown 
as contours in Fig.O As expected, the ability to detect a 
SSS is highly dependent on Nh and the source's effective 
temperature. Sources with low temperatures are more 
difficult to detect if significant Nn-absorption takes place. 
Without long (> fOksec) exposures, Nh ~ 10 22 cm~ 2 is 
the highest number that allows a detection of SSS emis- 
sion at 1 kpc distance with Swift. Since most novae are 
further away than 1 kpc, the situation is even worse. 

2.3. Targets 

Although primarily a mission to study the temporal 
evolution of 7-ray bursts (GRB), Swift also provides Tar- 
gets of Opportunity (ToO) observations for non-GRB 
events. Through ToO observations we have obtained X- 
ray observations of 9 recent (< f — 2 years since outburst) 
and 3 older (> 3 years) CNe. We consider all observa- 
tions, although those of V1188Sco on 21 May, 2006 and 
of V574Pup taken on 26 July, 2005 were extremely short 
with hardly any signal. Table Q] provides a list of impor- 
tant physical parameters associated with these novae in- 
cluding the maximum visual magnitudes, the date of dis- 
covery, the time to decline 2 magnitudes from maximum 
(t2), the early measured expansion velocities, E(_B — V) 
reddening values, Hi column densities, and distances. In 
the appendix lA.il we summarize background information 
for each nova in our sample. 

3. ANALYSIS OF THE Swift XRT DATA 

Currently the Swift archive contains 27 observations 
of 12 novae obtained with the X-ray Telescope (XRT) 
instrument (see Sect. 12.1] ). We analyzed the data using 
the Swift reduction packages in HEAsoft version 6.06 8 
and the latest calibration data (version 20060427). We 
started with the photon event files (level 2) that are avail- 
able in the public Swift archive and carried out three sta- 
tistical tests to determine if the sources were detected. 
We then calculated count rates or upper limits in cases 
where no detections can be claimed. The first criterion 
is the formal detection likelihood from comparison of the 
number of counts measured in the source extraction re- 
gion (circular with 10 pixels radius around expected sky 
position) with the background extraction region (annu- 
lus around source extraction region with inner radius 
10 pixels and outer radius of 100 pixels). We next de- 
termined if there was a concentration of counts near the 

8 http:/ /swift. gsfc.nasa.gov/docs/software/lheasoft/ 



center of the source extraction region. Finally, we com- 
pared the spectral distributions of photons in the source 
and background extraction regions, but this is only a soft 
criterion, as similar spectral distributions do not rule out 
the presence of a source. We present our st atist ical meth- 
ods in more detail in the appendix Sect. IA.2I to justify 
our tests. 

In Table [2] we provide all the Swift observational data 
including observation date, time since visual maximum, 
exposure time, number of counts in the source and back- 
ground extraction regions, and either the net source 
count rate per detect cell with ler uncertainty, or the 
upper limit. 

3.1. Non-detections, failure of all criteria 

Only upper limits could be established for the observa- 
tions of V476Sct, V1188Sco, LMC2005, and V1494Aql. 

In Fig. [3] we show the photon plot centered around the 
expected source position of V476 Set and LMC 2005. The 
extraction regions for source and background, as defined 
in Sects. l3l and IA.2I are marked with dark and light gray 
colors, respectively, and the counts within the source ex- 
traction region are plotted with small black bullets. In 
the upper right corner we list: the date of observation, 
source extraction region, the number of counts in the re- 
spective extraction regions, and the exposure times. In 
the upper left corner we display the source net count 
rate (or the 95.5-% upper limit if no significant excess 
count rate is found) as well as the detection likelihood, 
the number of net source- and background counts S and 
B, respectively, as defined in Eq. IA1I For V1188Sco, 
all criteria fail for all observations, and the source is 
clearly not detected. A similar case is encountered for 
all three observations of LMC 2005, where we show the 
last observation in the right panel of Fig. [3] Further, we 
found no significant detections in all three observations of 
V1494Aql. We found only a marginal detection in the 
Chandra observation of V1494 Aql five years prior to the 
Swift observations (see Sect. 13. 4| ). We used the proposal 
planning tool PIMMS in order to convert the Chandra 
count rate into an expected Swift XRT count rate. With 
various model assumptions PIMMS predicts more than 
0.04 XRT counts per second, while we measured less than 
0.0014 counts per second in each observation (Table [2]). 
This implies that the source has faded by at least a factor 
30 within five years. 

3.2. Weak sources < 0.01 cps 

VI 663 Aql was the first nova observed by Swift. Unfor- 
tunately, the first planned observation was truncated by 
a GRB after only 1.25 ksec. Nevertheless, we find a clear 
detection and a strong concentration of counts towards 
the center (Table [2]). The source spectrum is similar to 
the background spectrum, but that does not rule out a 
detection (see Fig. 2]). The spectrum is extremely weak 
and no conclusion can be drawn. A second observation 
seven months later was longer but yielded only an upper 
limit which implies a fading of the source. Neither of our 
other two criteria supported a detection. 

For V5116Sgr we only found a marginal (formal 93- 
percent) detection with a concentration of photons to- 
wards the center. The spectrum is too weak to apply 
our source-background comparison criterion (see Fig. [5]). 
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TABLE 1 
Basic nova parameters 



Object 


v 


Datc a 


t 2 b 


Vex C 


E(B-V) 




Dis. 


Rcfs. c 




(mag) 




(days) 


(km sec ) 


(mag) 


(cm -2 ) 


(kpc) 




V574 Pup 


8 


11/26/2004 


13 


650H,860P 


0.5 


6.5e21 


3.20 


1,2,3 


V382 Nor 


9 


03/19/2005 


13 


HOOP 


0.6-1.1 


1.8e22 


13.80 


3,3,3 


V1663 Aql 


10.5 


06/10/2005 


16 


700P 


~2 


1.8e22 


5.5±1 


4,5,4 


V5116Sgr 


<8 


07/04/2005 


20 


1300P 


0.24±0.08 


1.6e21 


11.30 


3,3,3 


V1188 Sco 


8.9 


07/26/2005 


12 


1730H,4000Z 


~1? 


5.0e21 


7.5 


3,3,3 


V1047 Cen 


8.83 


09/04/2005 


? 


850H 


3.3? 


1.6e22 


? 


-,3,- 


V476 Set 


10.9 


09/30/2005 


15 


1200H 


1.9±0.1 


l.le22 


4±1 


6,6,6 


V477 Set 


10.4 


11/10/2005 


3 


2700H,6000Z 


>1.3 


4.8e21 


11 


7,7,7 


LMC 2005 


12.6 


11/27/2005 


? 




0.15 


6.3e20 


50 


3,3,3 


V723 Cas 


7.1 


12/17/1995 


slow 


700H 


-0.5 


2.4e21 


4 


-,8,9 


V1494 Aql 


4 


12/03/1999 


6.6±0.5 


1300H,1850P 


0.6±0.1 


4.2e21 


1.6±0.2 


10,-, 11 


V4743 Sgr 


5 


09/20/2002 


9 


2400H 


low? 


1.4e21 


6.3 


3,-, 12 



a Date of visual maximum (mm/dd/yyyy). b Time to decline 2 magnitudes from maximum. c Expansion 
velocity. Taken from early IAU circulars. Trailing letters indicate how the velocity was measured: H 
= Full- Width at Half-Maximum, P = P Cygni absorption, Z = Full- Width at Zero-Intensity.* 1 The Hi 
column densities were obtained from the Heasarc NH tool. They are the line-of-sight column densities 
through the entire Galaxy toward the coordinates of each object. References for t2, E(B— V), and Dis- 
tance; 

1 = ISiviero et al.l 120051 ); 2 = R.J. Rudy 2006, private communica tion; 3 = t 2 : from AAVSO light curve 
and IAU circulars. E(B-V): using van den Bcrgh & Younger (1987) and the (B— V) color evolution or as- 
sumin g E(B -V) ~ Nff/4.8xl0 21 c m" 2 . Distance: using the My vs t 2 relationship o f Delia Valle & Livio 
1T995I); 4 =ILane et alJ 120061) ; 5 =IPuetter et al.l 12003) : 6 = IMunari et al.l l2006H)lPerrv et al.l 120051); 
7 =IMunari et al.l 12006al): 8 =lMunari et al.l 119961): 9 = llijima et al.1119981) ; 10 = IKiss fc Thomson! 
120001) ; 11 = IIiiima fc Esenoglul 120031) ; 12 = ILvke et aLl 120021) . 
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Fig. 3. — V476 Set (left) and LMC 2005 (right) were not detected. We show the recorded photon positions and mark the source extraction 
region with dark shading and the background extraction region (annulus with inner radius 10 pixels and outer radius 100 pixels) with light 
shading. The pixel- and sky coordinates are given in the bottom. The upper limits are given in the upper left. Observation date, extraction 
radius, number of source- and background counts per detect cell, and exposure time are given in the upper right. The percentage of source 
counts that reside within 25 percent area (radius 5 pixels) is given to the right of each source extraction region to give an impression of the 
concentration towards the center within the source extraction region. 
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Fig. 4. — Two Swift observations of V1663 Aql. The left panels (image) are the same as described in Fig. [3] The right panels provide the 
spectral information of the source (black) and background (downscaled to the area of the source extraction region; light shadings). The 
likelihood for the background spectrum being softer than the source spectrum is estimated by drawing 1000 random sample spectra with 
the number of source counts (in the upper panel 9) out of the background (0.31 X 99 = 31 counts) and counted the number of cases where 
the median energy of each sample spectrum was lower than that of the source spectrum. Much different values than 50 percent are a soft 
criteri on fo r a source detection. In the upper panel we find a clear detection and give the count rate with 1-a unce rtainty, AC according 
to Eq. IA6I and the number of net counts for source and background S and B (per detect cell) according to Eq. I All 
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TABLE 2 

Measured count rates in our sample 



Start 


At a 


Exp. 


Src. 


Bg. 


Rate b ' d 


Rl/2 C 


Date 


(d) 


(fa) 


(cts) b 


(cts) b 


(cts/ks) 


(%) 



V574 Pup 



05/20/2005 


175 


1.0 


10 


0.4 


12.0 ± 4.0 


60 


05/25/2005 


180 


1.9 


22 


0.1 


14.0 ± 3.0 


77 


05/26/2005 




0.01 














07/29/2005 


245 


1.1 


15 


0.3 


16.0 ± 4.0 


67 


07/30/2005 


246 


7.1 


51 


2.1 


8.5 ± 1.3 


73 


08/06/2005 


253 


7.8 


57 


2.0 


8.7 ± 1.2 


67 


08/09/2005 


256 


2.2 


22 


0.6 


12.0 ± 3.0 


64 


08/10/2005 


257 


2.8 


18 


0.8 


7.8 ± 2.0 


67 


08/11/2005 


258 


1.8 


11 


0.5 


7.3 ± 2.5 


56 


08/17/2005 


264 


4.7 


31 


1.2 


8.6 ± 1.6 


79 



V382 Nor 

01/26/2006 313 72 Zo 14.0 ± 2.0 75 



V1663 Aql 

08/15/2005 66 L3 9 03 8.6 ± 3.3 78 

03/04/2006 267 6.5 3 1.1 < 0.5 d 



V5116Sgr 

08/29/2005 56 O 5 2~0 1.2 ± 1.0 80 



VI 188 Sco 

05/21/2006 98 0.2 1 0.1 < 7 d 

06/17/2006 124 5.0 3 2.3 < 9.1 d 33 



V1047 Cen 



09/11/2005 66 3.9 19 1.3 5.6 ± 1.5 53 

01/23/2006 141 5.3 53 3.1 11.0 ± 2.0 72 



V476 Set 

02/12/2006 135 3.8 2 1.0 < 8.2 d 50 



V477 Set 



03/07/2006 117 6.2 15 1.8 2.7 ± 0.8 87 

03/15/2006 125 4.7 8 1.3 1.8 ± 0.8 50 



LMC2005 



12/02/2005 


5 


0.1 





0.1 


< 


1.0 d 





12/01/2005 


1 


1.1 


1 


1.3 


< 


8.6 d 


100 


12/02/2005 


5 


3.9 


2 


1.2 


< 


8.7 d 


100 


12/03/2005 


6 


5.7 


1 


1.7 


< 


7.0 d 


50 



V723 Cas 

01/31/2006 3698 6^9 147 L7 26.0 ± 2.0 69~ 
V1494 Aql 

03/10/2006 2289 1.4 0.3 < 1.4 d 

03/13/2006 2292 1.2 0.2 < 1.4 d 

05/19/2006 2359 2.5 1 0.5 < 1.0 d 100 

V4743 Sgr 

03/08/2006 1265 5^8 267 2~3 27.0 ± 2.0 70" 



a Interval between visual maximum and the XRT observation 
in days b per detect cell (10 pixels radius); count rates with la 
uncertainties fraction of counts within 5 pixels (%) d 95.5-% upper limit 



Two observations of V477 Set were carried out 8 days 
apart and resulted in two detections with a likelihood 
> 99 percent (see Fig. [SJ . The first observation suggests 
a weak source with a clear concentration towards the 
center and a much softer spectrum than that of the in- 
strumental background. The second observation reveals 
a lower count rate, lower detection likelihood, and a less 
clear concentration of photons towards the center. While 
this implies variability or that the source has faded over 
the eight day interval, within the 1-a errors, the count 
rate is consistent with that of the first observation. Nei- 
ther spectrum allows quantitative analysis. 

3.3. Strong sources > 0.01 cps 

V574Pup, V382Nor, V1047Cen, V723Cas, and 
V4 743 Sgr were detected by Swift (see Figs. I7HTT]). The 
spectrum of V382 Nor is hard whereas the background 
is softer. There appears to be some excess emission at 
wavelengths of prominent emission lines (e.g., 0.92 keV, 
Ne IX; 1.02 keV, Ne X; and 0.83 keV, Fe XVII) but a firm 
identification of these lines is not possible with the ex- 
isting data. There is also emission at energies with no 
well-known emission lines, particularly at high energies. 
The observed spectrum is too faint to derive any quan- 
titative conclusions. 

We obtained very clear detections of V1047Cen. In 
the first observation of November 2005 the photons are 
clearly concentrated towards the center, and the energy 
distribution is very different from the soft background. 
Two months later the detection is even stronger and the 
spatial as well as spectral distributions clearly support a 
detection. The count rate doubled from 5.6ctsksec _1 on 
9 November, 2005 to llctsksec^ 1 on 23 January, 2006. 
The spectrum appears hard, but there are too few counts 
to derive any spectral properties. 

The remaining novae were sufficiently bright to obtain 
spectral confirmation that they were detected during the 
SSS phase. More detailed analysis of these spectra, in- 
cluding additional Swift observations and supporting op- 
tical and near-IR spectra, is left for future papers. Here 
we only provide a brief analysis of the current Swift data 
sets. 

V574Pup was first observed in late May 2005. It was 
then observed multiple times from the end of July to mid 
August 2005 as it was perfectly positioned to cool the 
XRT detectors between GRB observations. In all nine 
observations, it was clearly detected with a high degree 
of significance (see Table [2]). V574Pup had an average 
of 7 — 14ctsksec _1 in each observation. The individual 
observations were too short for analysis, so we combined 
them, yielding a higher signal-to-noise ratio. We grouped 
the spectra by summing the number of counts and the 
number of background counts in each spectral bin. We 
calculated exposure time-weighted effective areas from 
the individual extracted effective areas. 

However, during the 50 days between the May 26 and 
July 29 observations, significant evolution probably oc- 
curred which we would miss by combining all spectra to- 
gether. In order to search for any evolution, we grouped 
the first two observations from May 2005 and the re- 
maining observations from July and August 2005. These 
are shown in Fig. 1121 The summed spectra of the two 
early observations shows a uniform distribution with en- 
ergy and may be an underexposed emission line spec- 
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Fig. 5.— Same as Fig.Hbut for V5116Sgr. 

trum. Unfortunately, we don't have any more data from 
this time period to increase the signal-to-noise ratio. The 
later observations sum up to what appears to be a contin- 
uous spectrum ranging from 0.3-0.5 keV (25 A to 45 A). 
This spectrum resembles that of a SSS. We carried out 
a blackbody fit using the most recent XRT response ma- 
trices (version 8) and assuming an emitting radius of 
7000 km. We found a satisfactory fit (Fig. ^ with the 
parameters T off = (293 ± 10) x 10 3 K (25.2 ± leV) and 
Nh = (2.5 ± 0.6) x 10 21 cm~ 2 (1-a uncertainty ranges). 
We fixed the radius because no unique solution could be 
found with a variable radius. The uncertainties include 
variations of the radius between 7000 and 9000 km. 

The XRT spectrum of V4743 Sgr is faint, and no spec- 
tral features can be identified. We also found no features 
in the RGS observation taken in September 2004 (see 
below). While the XRT count rates of V4 743 Sgr and 
V723Cas are similar, the spectra are different (Fig. [TT]) . 
V723Cas shows a clear peak occurring at 0.4 keV (31 A). 
iNess et all (j2006b| ) reported that the spectrum resembled 
a SSS, and they were able to fit a blackbody model. We 
repeated these fits using the updated XRT response ma- 
trices (version 8) and show the best-fit model in Fig. [TJJ 
We fixed the radius at 7000 km but included variations 
up to 9000 km in the calculation of the 1-a uncertainty 
ranges. We find T cff = (371 ± 15) x 10 3 K (32 ± 2eV) 
and N H = (3.3 ± 0.4) x 10 21 cm~ 2 . V723Cas is now, 
after 11 years, the oldest Galactic nova that is still ac- 
tive. Such a long duration of the SSS phase is unusual 
for a Classical Nova, and there is a possibility that the 
nova has transitioned into a SSS such as Cal83. Our 
observations of this interesting source are continuing. 

3.4. Supplemental observations 



In addition to the Swift observations, we extracted 
unpublished X-ray observations from the Chandra and 
XMM-Newton archives. A series of Chandra observations 
of V1494 Aql were carried out and we extracted the last 
observation (LETGS, ObsID 2681) taken 727 days after 
outburst (2001, Nov. 28). The source was only detected 
in the zeroth order on the HRC-S detector, but no dis- 
persed spectrum could be extracted. We found 56 counts 
on a background of 10, corresponding to an HRC count 
rate of 0.02 ± 0.003 cts sec -1 and a 14-tr detection. Unfor- 
tunately, the HRC detector does not provide the energy 
resolution to construct a spectrum, and no spectrum can 
be constructed from the dispersed photons as the back- 
ground is too high. 

We further analyzed XMM-Newton observations of 
LMC2005 taken 243 days after outburst (ObsID 
0311591201, 2006, July 18) and of V4743Sgr 1470 days 
after outburst (ObsID 0204690101, 2004, Sept. 30). 
LMC 2005 was again not detected. 

From MOS1 we extracted 0.11 ± 0.002 cts sec -1 for 
V4 743 Sgr. The source was also detected in the RGS 
with a count rate of 0.01 ± 0.0005ctssec _1 . Grat- 
ing spectra allow the calculation of a flux without the 
need of a model, and we measured a flux of 3 x 
10~ 13 erg cm- 2 sec" 1 from the RGS1 and RGS2 9 . Fig.[T5] 
shows a comparison of these last XMM-Newton spectra 
with the recent Swift observation. 

4. DISCUSSION 

The sample of X-ray observations of novae presented 
in this paper demonstrates that Swift significantly in- 
creased the number of Galactic novae observed in X- 
rays. We provide a large sample taken with the same 

9 0.35— 1.8 keV; corrected for the chip gaps 
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Fig. 6. — Same as Fig.[4]but for the two observations V477Sct. 
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Fig. 7. — Same as Fig.|4]but for six usable observations of V574Pup. 
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Fig. 8. — Continuation of Fig. [7] 
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Fig. 9. — Same as Fig. [4] but for V382Nor. The spectrum appears to be an emission line spectrum, but we note that it is oversampled 
and the lines cannot be as clearly identified. 
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Same as Fig. [4] but for the two observations of V1047 Cen. 
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Fig. 11.— Same as Fig. H but for V723 Cass (top) and V4743 Sgr 
different. 

instrument, which reduces problems from cross calibra- 
tion. Also, Swift is capable of carrying out monitor- 
ing observations with relatively little effort, the excep- 
tion being classical novae with large column densities. 
In Sect. 12.21 we demonstrated that one can not obtain 
enough counts in a typical ToO exposure ( 3-5 ks) when 
the column density exceeds N# > 10 22 cm~ 2 to observe 
an SSS phase. In this survey the estimated column densi- 
ties of four novae, V382Nor, V1047Cen, V1663Aql, and 
V476 Set, exceeded this value and only very hard spectra 
or non-detections were recorded. For Galactic novae with 



(bottom). While the count rates are similar, the spectra are very 

low column densities Swift is an excellent instrument for 
following all the stages of evolution in X-ra ys (e.g., the 
recurrent nova RS Qph lOsborne et aLll2006f ). 

Shortly after the outburst, X-ray emission from the 
WD is not expected to be observed. X-ray detections 
during this phase thus originate from other processes, 
e.g., shocks within the expanding shell or a shock set off 
by a collision between the expanding shell an d the atmo- 
spher e of the companion (as in RS Oph, e.g. JBode et al.1 
I2006TJ . Few o bservations during this phase have been car- 
ried out (e.g. iMukai fc~I shida 2 0011 ) and they only pro- 
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vide information on dynamics within the ejecta. With 
the flexibility of Swift a better assessment of pre-SSS X- 
ray emission can be achieved. 

The SSS phase is the brightest phase of X-ray emission 
(at constant bolometric luminosity), making it the best 
time to observe a nova. We found two clear detections of 
novae in this phase. With these two new SSS spectra, the 
number of Galactic novae clearly observed in this state 
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Fig. 15. — Comparison of the XRT spectrum of V4743 Sgr (top) 
with XMM-Newton RGS1 and RGS2 spectra from September 2004. 
All are consistent with a fading nebular spectrum, but the signal 
to noise and the resolution are insufficient to characterize the spec- 
trum as hard continuum or an emission line spectrum. 



is seven (|Orio et alJl2001bHDrake et alJl2003t iNess et al.l 
I2003L and reference therein). For V574Pup we were even 
able to identify the transition from an early emission line 
spectrum to a bright SSS spectrum within two months. 
V574Pup entered its SSS phase within 250 days after 
outburst. This is consistent with the evolution of V1974 
Cyg, V1494 Aql, and V4743Sgr (see Fig. [HI). The detec- 
tion of a SSS in V723 Cas suggests that nuclear burning is 
still going on more than 11 years after outburst. It is not 
known how much longer it will remain a SSS. More de- 
tailed discussions on V723 Cas will be presented by Ness 
et al. (in prep). 

We detected V4743 Sgr two years after XMM-Newton 
had already measured a very low level of X-ray emission 
(see Sect. 13.4] ). This shows that the post-turn-off evolu- 
tion can last a few years. V1494Aql was not detected, 
but five years prior to our observation it had been found 
to have marginal emission and no SSS spectrum. 

In order to compare the evolution of the novae detected 
in our sample we plot their distance-scaled X-ray bright- 
nesses as a function of individual time after outburst in 
Fig. HU In the bottom panel we show the X-ray light 
curves (in the same time units) of five more Galactic no- 
vae, but observed by different missions. The brightness is 
given as unabsorbed lu minosity as extracted from direct 
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measurements of flux es (iGreiner et al.l|2003bl : 
2001a; iBalman et all 119981 : INess et al.l 12003: 
2005K As a general rule, any X-ray detection of no- 
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vae earlier than 100 days after outburst arises from the 
hard, shock-generated phase that can either be emission 
lines or possibly bremsstrah lung emission (|Orio et al.l 
2001at iMukai fc Ishidal l200U) . Likewise, any X-ray de- 
tections more than 1000 days after outburst come from 
the nebular phase, with the three exceptions of GQ Mus, 
LMC 1995, and now V723Cas. These gen eral rules also 
hold for novae in M31 (see figure 3 in iPietsch et alJ 
I2005T ). 

Finally, we can use the parameters derived from the 
model fits to the SSS spectra of V574Pup and V723Cas 
to check the Swift counts that we predicted from our at- 
tenuated Cloudy model of an "average" classical nova 
ejection in Fig. [21 (Sect. I2.2|) . From our blackbody fits we 
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Fig. 16. — Top: XRT count rates of all sources in our sample, plotted against the day of observation after their outbursts. The count rates 
are corrected for distance squared which is uncertain. Multiple observations of the same targets are connected with dotted lines. Bottom: 
Observations from various other missions. The observed X-ray luminosities (not corrected for absorpt ion, except for V197 4 Cyg, which is 
instea d rescaled) of recent well-observed novae with SSS phases. Similar plots have been created by Pictsch ct al. (2005); Pictsch ct al. 
j2006h . 



obtained effective temperatures of ~ 3 x 10 5 K (25 eV) 
for both novae, while the derived column densities were 
2.7 and 1.7 x 10 21 cm" 2 for V574Pup and V723Cas, 
respectively. With these parameters we find from the 
grid of Cloudy models (Fig. [2]) the predicted counts in 
a 1-ksec observation at a distance of 1 kpc of ~ 10 for 
V574Pup and - 100 for V723Cas. Scaling these pre- 
dicted counts by the distances in Table [1] and the ex- 
posure times in Tabled] (t ea; p(ksec)/D(kpc) 2 ) gives 75 
counts for V723 Cas and 7 counts for the 30 July, 2005 
observation of V574Pup. These values are consistent 
with the detections given the uncertainties in the dis- 
tances and the use of a generalized nova model (e.g. a 
luminosity of 10 38 e rg s" 1 , an ejected mass ~ 10~ 4 M Q , 
iSchwarz et al.l 12007). For "normal" classical novae, cal- 
culations like those done for Fig. [3] can be used as a tool 
to plan X-ray observations by providing first-order esti- 



mates of expected X-ray brightness levels during the SSS 
phase. 
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APPENDIX 
A. APPENDIX 
A.l. Targets 

In this section we summarize some background information of our targets in addition to that given in Table [TJ 
collected mostly from unreferenced literature and private communications. 
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A. 1.1. V574Pup 

V574Pup was independently discovered by Tago a nd Sakurai (Nakano et al. 2004|). The evolution of its early light 



curve and spectral energy distribution is provided in ISiviero et alj ( 20051 ). The FWHM of Ha in the early spectrum 
was 650knis~ 1 but with P-Cygni absorption extending to 860 kms -1 . Based on the intrinsic colors of the nova both 
at maximum and at t2, the same authors derived an extremely low reddening value of E(B-V)~ 0.05. However, this 
value seems too low due to the low galactic latitude (b ~ 2°) and the distance (~ 3.5 kpc) of this nova. Optical and 
near-IR spectra obtained one year after out burst showed that V574Pup had entered a coronal phase with lines of 
[S VIII], [S IX], [Si VI], [Si VII], and [Ca VIII] (|Rudv et al.l l2005). Analysis of the ratio of the O I lines in this data set 
and those obtained later implies a higher reddening of E(B-V)= 1.27 (Rudy, private communication). 

A. 1.2. V382Nor 

V382Nor was discovered prior to visual maximum on 13.3 Mar, 2005 (jLiller Il2005ah . Examining the AAVSO light 
curve, visua l maximum p robably occurred on 19 March, 2005 at V^ 9. The B-V color obtained on 20 March was 
0.8 ± 0.11 ([LillerJ 2005a), which gives an E(B-V ) of 0.57 ± 0.17 based on the intrinsic color at maximum (B-V 
-0.23 ± 0.06; Ivan den Bergh fc Youngerlll987f). A nothe r determination of the reddening uses the intrinsic (B-V) 



color at t 2 , -0.02 ± 0.04 (|van d en Berg h fc Younger! 119871) . The (B-V) color at t 2 was ~ 1.1, which implies E(B- 



V)~ 1.1. This hi gh a reddening estima te is supported by the saturated Na I interstellar lines in the spectra taken 
at the same time (|Ederoclite et alJ l2005). The P-Cygni profiles seen in the early spectra imply an average expansion 
velocity of 1100 kms -1 . 

A. 1.3. V1663Aql 

V1 663 Aql was discovered on 9.2 June, 2005 during routine All Sky Automated Survey patrols (jPojmanski fc Oksanenl 
2005). Spectroscopy obtai ned one day after maxi mum showed an extremely red continuum implying significant red- 
dening toward the source (jDennefeld et al.l 12005). The emission lines had narrow P-Cygni profiles with expansion 
velocities of order 700 kms -1 . By 14 November V1663Aql had entered its nebu lar phase with [O III] in the optical 
and the coronal lines of [S VIII] , [S IX] , [Si VI] and [Si VII] present in the near- IR dPuetter et al. 2005). The expansion 
velocity, as measured from the FWHM of the emission lines, was 2000 kms -1 dPuetter et al.ll2005h The ratio of t he O I 



lines also confirmed the large extinction, with a value as large as E(B-V)= 2 ( Puetter et al.ll2005T ). lLane et~aT1 (|2006h 
reported an angular expansion of 0.2 mas/day at 2.2 /im with the Palomar Testbed Interferometer. This corresponds 
to a distance of 5.5 ± 1 kpc assuming expan sion velocities between 700 — lQOOkms -1 . If the expansion v elocities are 
larger , as measured bv lPuetter et all (|2005l ). then the distance must also be larger than that derived by lLane et al.l 
(f2006T) . 

A. 1.4. V5116Sgr 

On 4 July, 2005 Liller discovere d V5116Sgr already on the decline. His first observations implied that visual 
maximum was brighter than 8 mag (Liller 2005bJ). T he P-Cygni line profile of Ha suggested an expansion velocity of 
1300 km s -1 from spectra taken a day after discovery ([Gilmore fc Kilmartinl 120051 ). They also reported a B-V color of 
+0.47 ± 0.02 on 5 July. Using the intrinsic (B-V) color at maximum for V5116 Sgr implies an E(B-V)= 0.24 ± 0.08. 

A. 1.5. V1188Sco 

VI 188 Sco was d iscovered in late July 2005 (jPojmanski et al.1 l2005af ). The visual maximum was measured by 
ICooper et al.l (|2005l ). The only spectral information in the literature was obtained within days after visual maximum. 
Optical spectra were consistent with a typical CO type nova with br oad Ha, FWHM and Full Width at Ze ro Intensity 
(FWZI) of 1730 and 4000 kms -1 , and Fe II P-Cyg emission lines (|Naito et al.l 120051 : 1 Walter et al1l2005h . A 3 - 14 
micron spectrum taken one day a fter visual maxim um displayed a featureless continuum consistent with the Rayleigh- 
Jeans tail of the Planck function (ISitko et al.ll2005[). T he Na I D line was observed with three absorption components 
and an equivalent width of 0.5 nm fWalter et al.ll2005T ) implying a large extinction. 

A. 1.6. V1047Cen 

Little is known about V1047Cen due to its extreme southern declination (> —62°). It was discovered bv iLillerl 
(|2005cD on 1.031 September, 2005. Visual maximum probably occurred three days later at V= 8.83 but there is very 
little reported data on the early light curve, so its subsequent behavior, including its t 2 decay time, is unknow n. The 
spectr um at maximum resembled that of V5114Sgr and V5116Sgr with an Ha FWHM of 840 ± 50 kms -1 (jLiller I 
120053) . 

A.1.7. V476Sct 

V476 Set was discovered after visual maximum on 30.5 Sept, 2005 by A. Takao (|Soma et al.ll2005h . Haseda (2005) 
reported an independent discovery on the same day with a visual magnitude of 10.9. Early optical spectra showed 
many emission lines, notably of Fe II (jMunari et al.ll2006bl ). In the spectra the prominent lines showed double-peak 
profiles with velocity separations of ~ 700 km s -1 and FWHM of order 1000 kms -1 . Optical and near-IR spectra 
obtained six weeks after maximum showed similar emission features but with a strong red continuum implying the 
presenc e of a dust shell in the ejecta. The reddening inferred fro m the ratio of the O l lines at that time was E(B- 
V)~ 2 ()Perrv et al.ll2005l ). This value is similar to that derived bv lMunari et al.l (|2006bl ) using the color evolution and 
strength of the diffuse 6614A interstellar band. 
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A.1.8. V477Sct 

The second nova discovered in the constellation Scutum in 2005, V477Sct, was detected independently by G. Po- 
jmanski and K. Haseda on images obtained on the 11th and 13th of O ctober, 2005, respective ly, and the maximum 
visual magnitude was recorded on 13.07 October at V = 10.44mag (|Poimanski et alj|2005b[ ). Early JHK spec- 
troscopy of V477 Set sh owed broad H I emission lines with FWZI of 6000kms _1 and no evidence of dust emission 
( Poimanski et al.1|2005bf) . Optical spectra confi rmed broad emission lines with a FWHM value of Ha = 2900 km s -1 
( Fuiii fc Yamaokal f2005) . The optical spectra of iMunari et al.l (|2006af ) showed that unlike V476Sct which was rich in 
Fe II lines, V477 Sct had He and N lin es. Although sparse, the light curve revealed a rapid decline with an estimated 
t 2 of only 3 days (jMunari et al.ll2006aD . This rapid decline was consistent with the large expansion velocities. Optical 
and near-IR spectra take n on 15.094 Novem ber, 2005 showed a continuum increasing toward the red, possibly due 
to thermal dust emission (|Mazuk et al.| [2005). Numerous emission lines were observed which had developed double 
peak profiles with a FWHM of 2700 km s -1 . The features present included the coronal lines [S VII] and [S IX] but not 
[Ca VII] or [Si VI]. The ratio of the O I li nes implied a substant ial reddening of E(B-V)= 1.2 (jMazuk et al.ll2005l ). A 
similar value of > 1.3 was determined by IMunari et a l. (2006a]) based on the early color evolution, extinction maps 
along the line of sight, and the large equivalent width of the diffuse interstellar band at 6614 A. 

A. 1.9. LMC2005 

Liller discovered LMC2005 on 26.16 November 2005, but it was als o present on images that he had taken four 
days earlier. Visual maximum occurred 27 November at V — 12.6 mag (iLiller Il2005dh . The light curve evolution was 
extremely slow and in early 2006 exhibited photometric evidence of dust formation (F. Walter, private communication). 

A.1.10. V723 Cas 

V723Cas is a slow nova. It was discovered on 24 August, 1995 by Yamamoto (Hirosawa etaTlll995h at V = 9.2 mag. 
Prediscovery observations showed that it had slowly brightened over the previous 20 days. The subsequent light curve 
was unusual. V723 Cas reached a maximum magnitude of V ~ 8.6 mag one hundred days after discovery. At that point 
V723Cas exhibited a flare to 7.1 mag on 17 December, 1995. D uring the flare the U — B colors became significantly 
bluer while the B-V color remained constant (Mu nari et al.|[l996T) . The light curve also had ad ditional secondary flares 
durin g the follo wing 400 days. V723 Cas has been extensively observed from the ultraviolet (|Gonzalez-Riestra et al.1 
1996) to radio (|Hevwood et al.|[2005f) wavelengths. Recent optical and NIR spectra showed that V723Cas reached 
a "coronal" phase with em ission from high ionization stages such as [Fe X] (6373 A) similar to spectra obtained of 
GQ Mus in the late 1980s dKrautter fe Wilriamslll989l ). The spectral similarity suggests ongoing nuclear burning on 
the WD (|Krautter fc Williams! 1 19891) . As a result, we requested Swift observations in order to search for the X-ray 
emission that would be related to ongoing nuclear burning. 

A. 1.11. V1494Aql 

V1494 Aql was discovered on 1.8 December, 1999 at V — 6.0mag bv lPereiral (|1999f ). Within a few days of discovery 
this nova reached a peak magnitud e exceeding 4. It subse quently declined by two magnitudes in 6.6 ± 0.5 days thus 
classifying V1494 Aql as a fast nova (|Kiss fc Thoms on 2000). Three Chandra ACIST observations were c arried spectrum 
with no obvious soft component. However, a bright SSS spectrum was present in the August spectrum (|Starrfield et al.l 
2000). One month late r, two grating obs ervations were obtained on 28 September and 1 O ctober, 2000, which showed 
remarkable variability (|Drake et al.l l2003). The spectrum resembles that of the SSS Cal83 (|Paerels et al.ll200ll ). except 
that it appeared somewhat hotter and the "emission features" were at different wavelengths. A final observation 
with Chandra LETGS one year later (ObsID 2681) showed no emission in the dispersed spectrum but a reasonable 
detection in the zeroth order on the HRC-S detector. We extracted 56 counts on a background of 10, corresponding 
to a count rate of 2cts/ksec and a 14-er detection likelihood. Unfortunately, the HRC detector does not provide the 
energy resolution to construct a spectrum. 

A.1.12. V4743Sgr 

V4743 Sgr was discovered bv lHaseda et al.1 (|2002l ). A distance of 6.3 kpc is assumed b ased on infrared observations 
(jLvke et al.l[20 02). Chandra grating observations were carried out on 19 March, 2003 (|Ness et~all[200l . The light 
curve showed large amplitude oscillations (~ 20% of the count rate) with a period of 1324 sec (almost 40% power) 
with the 2nd and 3rd harmonic overtones present. V4743Sgr was observed by Chandra on five more dates and details 
of the spectral evolution will be presented by Ness et al. (in prep). 

A. 2. Extraction of source counts 

In this section we give a detailed account of our methods to determine the count rates and detection probabilities. 
We determined the positions of the sources on the chip from the individual sky coordinates and used circular extraction 
regions (radius 10pixels= 23.6", which encircles ~ 80.5 percent of the total PSF). For the extraction of the background 
we defined annular extraction regions around the source position with an inner radius of 10 pixels and an outer radius 
of 100 pixels. We carefully checked to make sure that no sources were present in the background extraction regions. We 
scaled the number of extracted background counts by the ratio of extraction regions (99.0) to determine the number 
of expected background counts per detect cell. We assume Poissonian noise and thus apply a maximum likelihood 
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estimation in order to determine the net source count rate. We define two model parameters S and B that represent 
the number of expected counts of source and background for the source extraction region and calculate the number of 
expected counts in the source- and background extraction regions N s and Ni,: 

N s = aS + B N b = (1 - a)S + 0B (Al) 

with a = 0.80 the fraction of the PSF included in the detect cell and (3 — 99 the ratio of background- to source extraction 
areas. Under the assumption of Poissonian statistics we can calculate the probability of finding the measured numbers 
of counts n s and n b in the respective extraction regions when S and B are given: 

p= ^s e -N 3 ^ e -N b (A2) 

n s \ n b \ 

and thus the likelihood can be calculated 

C = -2 lnP = -2n s \n(aS + B) - (aS + B) (A3) 
+n b ln((l - a)S + (3B) - ((1 - a)S + @B) + const . 

We seek solutions for S and B where C reaches a minimum, which holds for N s — n s and N b = n b , and thus 

n b ~(3n s n s -a(n s +n b ) 

b = i a' B = — i a — ^ A4 ) 

1 — a — ap 1 — a — ap 

leading to a minimum value C m in of 

£ m i„ = -2[n s (hxn s - 1) + n b (lnn h - 1)] . (A5) 

With the definition of £ the range between C m i n and C m i n + 1 is equivalent to the 68.3-percent uncertainty range of 
the critical parameter S. We thus obtained the formal l-a errors by varying S until the value of C had increased by 
1.0 from the minimum while leaving B fixed at the value obtained from Eq. IA41 

The value of £ m i n can also be used to calculate the detection probability (in percent) based on the likelihood ratio 
test, which is a statistical test of the goodness-of-fit between two models. A relatively more complex model is compared 
to a simpler model (null hypothesis) to see if it fits a particular dataset significantly better. We define the likelihood 
for the null hypothesis Co as C(S — 0), and 

AC = C(S) - C (A6) 

quantifies t he improvement in C after including S > 0. The parameter B is an "uninteresting parameter" as defined by 
lAvnil (|1976f ). and we thus convert AC to the probability for one degree of freedom using the IDL 10 function chisqr_pdf . 

In cases where no counts were found in the detect cell (n s = 0), or when the number of counts in the detect cell is 
smaller than that expected from the background (n s < n b ), we calculated the 95-percent upper limits from the number 
of expected background counts and those actually measured in the detect cell. We again used the likelihood ratio test 
and solved Eq. IA6I for C(S) according to Eq. IA3I yielding the value of S that returns AC — 4. We also calculated 
upper limits in cases where the derived l-a uncertainties were larger than the measured count rate. 

In addition to the formal detection likelihood we consider whether the photons inside the source extraction region 
show some concentration towards the center, and whether their energy distribution is different from that of the 
instrumental background. In order to assess the concentration towards the center, we give the percentage of counts 
that we found in the central quarter of the detect cell in the last column of Table El From the PSF we expect 
60.5 percent of source photons to be within a circle of radius 5 pixels, and any number much below this fraction is 
likely not a source. We also studied the spectral energy distribution of the background compared to that of the source 
as an auxiliary criterion. We computed the median of recorded energies of the source counts and compared with 
the same value from the background. These values are marked by vertical lines (gray dotted for background and 
black dashed for source in Figs. 131 to [TTj) . In order to assess the probability that the background spectrum is softer 
than the source spectrum, we generated 1000 spectra with n s counts drawn as random sub-samples out of the pool 
of extracted background counts n b and calculated the median value for each generated spectrum. While a situation 
in which 50 percent of all random cases returned lower median energies does not imply a non-detection, any strong 
deviation from 50 percent is supportive of an independent source spectrum. 
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